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We have recently identified metal-sandwich (MS) crystal structures and shown with ab initio 
calculations that the MS lithium monoboride phases are favored over the known stoichiometric ones 
under hydrostatic pressure [Phys. Rev. B 73, 180501 (R) (2006)]. According to previous studies 
synthesized lithium monoboride tends to be boron-deficient, however the mechanism leading to this 
phenomenon is not fully understood. We propose a simple model that explains the experimentally 
observed off-stoichiometry and show that compared to such boron-deficient phases the MS-LiB 
compounds still have lower formation enthalpy under high pressures. We also investigate stability 
of MS phases for a large class of metal borides. Our ab initio results suggest that MS noble metal 
borides are less unstable than the corresponding AlB2-type phases but not stable enough to form 
under equilibrium conditions. 
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1. INTRODUCTION 

The interest in the AIB2 family of metal diborides 
re-emerged after the discovery of superconductivity in 
MgB2 with a surprisingly high transition temperature 
of 39 K0. Boron p-states have been shown to be 
key for both stability and superconductivity in these 
compounds 0, 0, 01- MgB 2 is a unique metal diboride 
because it has a significant density of boron pa-states 
at the Fermi level which give rise to the high T c super- 
conductivity, and yet enough of them are filled for the 
compound to be structurally stable 0,0)0- The effec- 
tively hole-doped noble- and alkali-metal diborides would 
have higher pa density of states (DOS) at Ef, but they 
have been demonstrated to be unstable under normal 
conditions 0| . The effort to achieve higher T c has thus 
primarily focused on doping magnesium diboride with 
various metals: however, doping this material has proven 
to be difficult [5J] and no improvement on T c has yet been 
reported. According to a recent theoretical study of non- 
local screening effects in metals, MgB2 may already be 
optimally doped0]. Lithium borocarbide with a doubled 
A1B 2 unit cell has been suggested as a possible high-T c 
superconductor under hole-doping 7] , but disorder in the 
heavily doped Li^BC appears to forbid superconductivity 
above 2K@. 

In this work we investigate whether there could be sta- 
ble high-T c superconducting metal borides in configu- 
rations beyond the standard AIB2 prototype. We have 
recently proposed metal-sandwich (MS) structures MSI 
and MS2, which also have sp 2 layers of boron but sepa- 
rated by two metal layers 9] . Despite their rather simple 
unit cells these structures have apparently never been 
considered before. As we demonstrate below, identifi- 
cation of the MS structures is not straightforward be- 
cause they represent a local minimum not usually ex- 
plored with current compound prediction strategies 
We reveal trends in the cohesion of MS phases by calcu- 
lating formation energies for a large class of metal borides 
and show that some monovalent-metal borides benefit 



from having additional layers of metal. The MS noble- 
metal borides still have positive formation energy, but 
they are less unstable than the AlB2-type phases. This 
result helps resolve the question of what phases would 
form first in the noble-metal boride systems under non- 
equilibrium conditionspiJ El El El 



Our main finding concerns the Li-B system, in which 
the MS lithium monoboride is stable enough to compete 
against the known stoichiometric phases. According to 
our previous ab initio calculations the MS lithium mono- 
boride is comparable in energy to these phases under nor- 
mal conditions, but it becomes the ground state at 50% 
concentration under moderate hydrostatic pressures 0. 
Here we extend the analysis to non-stoichiometric Li-B 
phases which could potentially intervene in the synthe- 
sis of the MS phases. In particular, synthesized lithium 
monoboride with linear chains of boron is known to be 
boron-deficient for reasons not fully understood so far. 
We simulate the incommensurate LiBj, phases (notation 
explained in Ref. El) by constructing a series of small 
periodic Li2„B m structures and show that the minimum 
formation energy is achieved for y « 0.9, in very good 
agreement with the experimentally observed values. Us- 
ing this simple model of the off-stoichiometry phases with 
linear chains of boron we demonstrate that relative to 
them MS-LiB still has lower formation enthalpy under 
high pressures. Simulations of other alkali-metal borides, 
MB y (M = Na, K, Rb, Cs), suggest that that these nearly 
stoichiometric phases might form under moderate pres- 



The paper is divided in the sections describing: 2) sim- 
ulation details; 3) construction of the MS prototypes; 4) 
stability of MS phases for a large class of metal borides; 5) 
detailed investigation of the Li-B system; 6) simulations 
of other monovalent and higher-valent metal borides; 7) 
summary of the electronic and structural properties of 
the MS phases. 
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2. COMPUTATION DETAILS 



Present ab initio calculations are performed with Vi- 
enna Ab-initio Simulation Package VASP 0, 0] with 
Projector Augmented Waves (PAW) ^3] and exchange- 
correlation functionals as parametrized by Perdew, 
Burke, and Ernzerhof (PBE) 18] for the Generalized Gra- 
dient Approximation (GGA). Because of a significant 
charge transfer between metal and boron in most struc- 
tures considered we use PAW pseudopotentials in which 
semi-core states are treated as valence. This is espe- 
cially important for the Li-B system as discussed in Refs. 
0, Il9| . Simulations are carried out at zero tempera- 
ture and without zero-point motion; spin polarization is 
used only for magnetic alloys. We use an energy cut- 
off of 398 eV and at least 8000/ (number of atoms in unit 
cell) k-points distributed on a Monkhorst-Pack mesh ^2(^ . 
We also employ an augmented plane-wave-|- local orbitals 
(APW+lo) code WIEN2K to plot characters of electronic 
bands |2l|. All structures are fully relaxed. Our care- 
ful tests show that the relative energies are numerically 
converged to within 1~2 meV/atom. 

Construction of binary phase diagrams A^Bi-^ is 
based on the calculated formation enthalpy Hf, which 
is determined with respect to the most stable structures 
of pure elements. For boron there are two competing 
phases a-B and /3-B|22j; we use a-B (Ref. 4j), theoreti- 
cally shown to be the more stable phase at low temper- 
atures and high pressures A structure at a given 
composition x is considered stable (at zero temperature 
and without zero-point motion) if it has the lowest for- 
mation enthalpy for any structure at this composition 
and if on the binary phase diagram Hf(x) it lies below 
a tie-line connecting the two stable structures closest in 
composition to x on each side. 



3. IDENTIFICATION OF MS PROTOTYPES 

Data-mining of quantum calculations (DMQC), intro- 
duced in our previous workj2^|, is a theoretical method to 
predict the structure of materials through efficient re-use 
of ab initio results. The DMQC iteratively determines 
correlations in the calculated energies on a chosen library 
of binary alloys and structure types. The last work has 
demonstrated that for a set of 114 crystal structures and 
55 binary metallic alloys the method gives an almost per- 
fect prediction of the ground states (within the library) 
in a fraction of all possible computations [U Eif . The 
speed-up (commonly by a factor from 3 to 4) is achieved 
by the method's rational strategy for suggesting the next 
phase to be evaluated. An essential feature of these calcu- 
lations is the full relaxation of the considered structures, 
which ensures an accurate determination of the correla- 
tions in the chosen library 0. 

We have recently begun expanding the 114 x 55 library 
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FIG. 1: (color online). Energy barrier as a function of 
boron fractional coordinate ds for a transition path from V2 
(d B « 0.25 26|) to MSI (d B » 0, see text) for different metal 
borides. For each value of ds the energy is minimized by ad- 
justment of the remaining free parameters o, c, and the metal 
fractional coordinate cLm- 



of ab initio energies of binary alloys [23j to include metal 
borides. Boron tends to form covalent bonds in inter- 
metallic compounds; to have this correlation in future 
predictions with the DMQC we needed first to add a few 
compound- forming metal-boride systems into the library. 
Introduction of a new system involves calculations of en- 
ergies for all the prototype entries in the library. Sur- 
prisingly, in the very first system considered, Mg-B, one 
of the fee structures with 4-atom unit cell at 50% con- 
centration, A2B2 fee- (111) (or V2 Hi|), relaxed almost 
all the way down to the AlB2-MgB2^hcp-Mg tie-line. 
Significant relaxations are not uncommon in our simu- 
lations; they usually correspond to the transformation 
from a starting configuration to a known stable proto- 
type and are automatically detected by the change in the 
symmetry. The magnesium monoboride phase, however, 
retained its original space group R3m (#166). 

Having examined the relaxation process we found that 
there is a continuous symmetry-conserving path from V2 
to a new structure MSI 0, E3| • V2 has 4 atoms per unit 
cell with 4 free parameters a, c, ds 1 and oIm (the last 
two are fractional distances between boron and metal 
layers), so that atoms are constrained only to the verti- 
cal lines (see Fig. 1). The perfect fee lattice corresponds 
to c/a=4-^/2/3 and oIb = d^/=0.25, but this special case 
does not grant additional symmetry operations and local 
relaxation have been seen in some metallic systems |26j. 



In metal borides a more dramatic transformation leads 
to a much more stable configuration: boron atoms re- 
arrange themselves to form covalent bonds in a hexago- 
nal layer (alg — * 0) rather then share electrons in close- 
packed triangular layers, while metal atoms remain in a 
close-packed bilayer. We have checked other alkali, al- 
kaline and transition metal borides not present in the 
DMQC library and confirmed that they all benefit from 
this transformation; however, some electron-rich systems 
might not escape from the local fcc-type minimum, as 
shown in Fig. 1 for FeB. This could be a reason why 
the MSI prototype has apparently been overlooked so 
far. We would like to point out that identification of new 
prototypes is not an intended function of the DMQC. 
This interesting accidental result should be credited to 
the exhaustive consideration of all candidates (regardless 
of how unlikely they seem to be a stable phase - an fee 
supercell is hardly a suitable configuration for a magne- 
sium boride phase) and the careful structural relaxation 
in the calculation of their ground state energies. 

We proceed by constructing a library of related MS 
prototypes. Structures, where the metal atoms closest 
to boron sit directly above the center of boron hexagons, 
are uniquely specified by the positions of the metal lay- 
ers (such as a, (3, 7, in Fig. 2). The MSI structure can 
thus also be labeled as \a(3\: the Greek letters show the 
positions of the two metal layers and vertical bars cor- 
respond to boron layers. A hexagonal supercell for this 
phase |27j is obtained when the last metal layer matches 
the first: |a/3|/37|7a|. The fourth metal layer can alter- 
natively be shifted back to site a (see Fig. 2(d)), result- 
ing in another structure at the same stoichiometry MS2 
(|a/3|/?a|)H|2|. Examples of more metal-rich structures 
are \a(3a\ (MS3)[ii| and |a/?7|. Various stoichiometries 
can also be achieved by combination of smaller cells, 
i.e. \aP\/3\ (MS4) and |a|a/3|/3a|. In this notation the 
AlB2-prototype is labeled simply as \a\, which we will 
use henceforth to avoid confusion with the aluminum di- 
boride compound. Positioning metal atoms above boron 
hexagon centers is not the only possibility. Stac king 
faults have been experimentally observed in MgBo[3f|. 
though such defects have been shown to be energetically 
costly for this compoundlU- We have constructed a pe- 
riodic structure (S-MB2) with 3 atoms per unit cell where 
the metal atoms in |a| are shifted along (a+b) to be 
above the middle of a boron-boron bondl 



4. STABILITY OF THE MS PHASES 

Formation of a particular compound in systems with 
a few competing phases is determined by a number of 
factors, i.e. the ground state energy, synthesis conditions, 
thermodynamic and kinetic effects. Comparison of high- 
throughput ab initio results[3ll3jl3 with experimental 
databases has shown that the calculated total energies 
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FIG. 2: (color online). Known AIB2 and proposed MS struc- 
tures for metal borides. The hexagonal layers of boron (grey) 
are separated by triangular layers of metal (yellow, only four 
atoms per layer are shown). Except for the 8 structure shown 
in b') |32|l . all structures have metal layers positioned above 
the middle of the nearest boron hexagons. The alternative 
notations given in brackets are explained in the text. The 
concentration of metal in these structures is 33% (AIB2), 33% 
(S), 50% (MSI), 50% (MS2), 60% (MS3), and 43% (MS4). 



alone allow to identify the correct phases observed in the 
experiment in about 96.7% of investigated cases (Eq. 3 
in Ref. Hj|)- In this section we use the total energy 
criterium to narrow down the set of systems in which the 
MS phases might occur. 

We first calculate the formation energy for a large li- 
brary of alkali, alkaline, and transition metals in the \a\ 
and MSI configurations (Fig. |3J. Our results for |a|-MB2 
are consistent with the previous calculation by Oguchi[4J 
(note that we put noble metals in the first valence group) . 
The MSI phases exhibit a similar trend in cohesion: they 
are most stable for tetra-valent metal borides. It is con- 
venient to analyze the stability of the MS phases by com- 
paring them to the corresponding |a|-phases because the 
MS structures are effectively a combination of the |a|- 
structure and additional layers of metal. Three imme- 
diate effects can be expected from the insertion of an 
extra metal layer: 1) different strain conditions between 
the boron network and the triangular layers of metal; 2) 
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FIG. 3: (color online). Calculated formation energy for |q| 
(top) and MSI (bottom) metal borides. The results for \a\- 
MB2 are consistent with the previous study[J] (mind the mi- 
nus sign difference between formation energy and heat of for- 
mation) . 



different doping level of the boron layer; 3) significant re- 
duction of interlayer overlaps between p-orbitals of boron 
due to the increase in the interlayer distance. To decou- 
ple these effects we calculate relative stability of the MSI 
phase with respect to phase separation into \a\ and pure 
element for the large library of metal borides. The rel- 
ative stability for compound M^Bi^ (= X) is defined 
as 



AEx = Ex — — 



(l-x)E ]a{ 



(1) 



(all energies are per atom) , and reflects whether a metal 
layer prefers to be in a layered boride environment or 
stay in pure bulk structure. To illustrate the amount 
of strain in the system we plot this energy difference 
versus equilibrium intralayer distance in pure fee or hep 
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FIG. 4: (color online). Relative stability of the MSI metal 
borides with respect to phase separation into \a\ and pure 
metal phases (Eq. 1) as a function of the close-packed dis- 
tance defined in the text. The metal borides are grouped and 
labeled according to the valence of metal (from 1 to 4). Inset: 
the same for noble- and alkali-metal borides. 



bulk metal structure, whichever is more stable at zero 
temperature [s^j . Figure 4 shows that monoborides of 
metals in the same valence group and similar dimensions 
(for example Zn, Cd, and Mg) have close relative stabil- 
ity. Metal layers mismatched with the boron layer (for 
example Be, Ca, Na, and K) cause a significant energy 
penalty when inserted in the respective diborides. 

Another general trend captured in Fig. 4 is a con- 
sistent decrease of the MSI phase relative stability with 
the increase of valence electrons (up to four). Relative 
stability depends on the binding mechanisms in all the 
three phases in Eq. 1 and the analysis of its variation 
with the metal valence is not straightforward. For mono- 
valent metal borides the observed gain in binding for the 
MSI phase is consistent with the fact that diborides of 
low-valent metals have available pa bonding states and 
stabilize as the metal valence increases^ 0|- However, 
charge redistribution in these phases may follow different 
scenarios: in \a\ metal atoms are exposed to boron and 
become almost entirely ionized^, while in MSI boron 
extracts charge through the surface of the metal bilayer 
and likely leaves more charge in the metal system (Sec- 
tion 7). With increasing valence in the transition metal 
series of diborides bonding pa and d-p~n bands become 
occupied, so that the binding reaches its maximum for 
Group IV metals and eventually goes downQ. Conse- 
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quently, we observe a noticeable increase in the relative 
stability of the MSI structure for higher- valent metals (in 
fact, all metal diborides with at least five valence elec- 
trons benefit from insertion of an extra metal layer with 
the largest gain of -340 meV/atom obtained for RhB). 
However, these electron-rich systems allow other phases 
with significantly lower energies (prototypes NiAs, NaCl, 
FeB-b, etc.|l3)- Hence, we focus on low-electron systems 
that have been shown to stabilize through incorporation 
of extra metal layers and could compete with existing 
phases. 



5. LI-B SYSTEM 

Overview. A few compounds at different stoichiome- 
tries have been reported for the Li-B svstempli. [3(| I37I 
EM [hj [TIJ . On the boron-rich side the experimentally 
reproducible compounds LiaBi4 and LiBs have large unit 
cells with fractional occupancies ^3,HSH3 an d cannot be 
presently simulated with ab initio methods with desired 
degree of accuracy. The composition of the most lithium- 
rich LiBy compounds (near 50% concentration) appar- 
ently depends on synthesis conditions and post-synthesis 
treatment, as the reported values for y range from 0.8 to 
1 (notation explained in Ref. [13)- I n the early experi- 
ments the formed compounds were ascribed compositions 
Li 5 B 4 |43,|4^| or Li 7 B 6 [44|]; Wang at al. used a rhombohe- 
dral model to explain the observed x-ray patterns |42l |45| . 
However, a more consistent interpretation of the avail- 
able x-ray data on nearly stoichiometric lithium mono- 
boride has been recently given by Liu et aZ. 38]. The 
authors demonstrated that the main x-ray peaks can be 
indexed with a four-atom hexagonal unit cell a-LiB (Fig. 
5(a)), which consists of linear chains of boron embedded 
in hexagonal lithium shells 38] . 

While the simple a-LiB sheds light on what the struc- 
ture of the lithium monoboride is, an important ques- 
tion remains open as to why the LiBj, compounds are 
boron-deficient. Worle and Nesper have offered an 
insightful model of LiBj,, in which the boron chains 
are uncorrelated and incommensurate with the lithium 
sublattice^a]- By using a large unit cell containing 
32,000 atoms the authors reproduced a kink at 26 w 
60° in the x-ray pattern and attributed it to the av- 
erage boron-boron distance of 1.59 A. They also sug- 
gested that the boron chains could be dimerized or have 
vacancies 0. According to a recent theoretical study, 
boron chains in lithium monoboride are not expected to 
dimerize but might indeed be able to slide freely along 
the lithium sublattice|39|. 

Model of the LiB y compounds. Simulation of the dis- 
ordered LiB y compounds is essential for finding a pos- 
sible stability region of the MS-LiB phases. While the 
large unit cells with thousands of atoms are needed to 
reproduce the x-ray data, such sizes prohibit the use of 
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FIG. 5: (color online) . Structures with linearly chained boron 
immersed in a lithium sublattice: the known a-LiB and /3-LiB 
have commensurate sublattices; the proposed Li2nB 77i series 
(only LigB7 is shown) allows to simulate the incommensurate 
off-stoichiometry LiB y phases. 



ab initio methods for ground state energy calculations. 
Therefore, we simulate the incommensurate LiBj, com- 
pounds by constructing a series of relatively small com- 
mensurate Li2 n B m phases. The number of lithium atoms 
in a unit cell must be even since they occupy alternating 
sites along the c-axis. To determine the optimal rela- 
tive position of the two sublattices we fix only the z- 
components of one lithium and one boron atoms and al- 
low all the other degrees of freedom to relax. We find 
that the relative placement matters only for the smallest 
Li2B2 unit cells: as we have shown in Ref. [9| the energy 
difference between a-LiB and /3-LiB is 10 meV/atom. 
For all other periodic structures the barriers to sliding 
for the two sublattices are below 1 meV per unit cell. 
The situation is similar to the relative motion in multi- 
walled carbon nanotubes, where the rigid layers interact 
weakly with one another: in long-period commensurate 
nanotubes the barriers to intertube sliding are extremely 
small, and in incommensurate ones the intertube sliding 
mode is gapless|47|. Local relaxations in Li2 n B m (to > 2, 
to =/= 2n) are insignificant due to the rigidity of the boron- 
boron chains and a smooth charge density distribution 
along the chains |39j. 

Stability and structure of the LiB y compounds. Figure 
6 shows an immediate benefit for the lithium monoboride 
to change composition: as the level of lithium concen- 
tration increases by a few percent the phase undergoes 
stabilization by over 20 meV/atom at zero pressure. The 
points on the enthalpy versus concentration plot for the 
Li2«B m series nicely fit to a parabola (with minimum at 
y « 0.894 at zero pressure). This leads to an interesting 
situation, in which the LiBj, phases are actually stable 
over a range of concentrations. The lower boron concen- 
tration limit corresponds to y m in ~ 0.874 (the tangent 
to the parabola going through x = 1, shown in Fig. 6), 
while the higher one depends on the location of LiB3 on 
the phase diagram and should be around y max ~ 0.9. 
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FIG. 6: (color online). Calculated phase diagrams for the 
Li-B system. Hollow red points are M2 n B m phases, from left 
to right: a-LiB, /3-LiB, LiioBg, LisB7, LieB 5 , LiioBs, Li4B3 
for all pressures, LigB 5 for P = 12 GPa, and LigB 5 , Li2B for 
P = 30 GPa. The parameters of the parabolic fit are given 
on each panel; the solid lines are tangent to the parabolas 
at x = 0.534, x = 0.560, and x = 0.584 for 0, 12, and 30 
GPa pressure, respectively. |a|-LiB is the AIB2 prototype and 
B32 is the pseudodiamond struct ure|3S)|. Note the different 
enthalpy scale for different pressures. 



The allowed concentrations are in excellent agreement 
with the Worle and Nesper's value of y = 0.9 inferred 
from the analysis of the x-ray data^(| . Considering that 
the lithium and boron sublattices are nearly independent, 
it seems possible to manipulate the stoichiometry with 
active solutions by removing the alkali metal through 
the surface of the sample. By using tetrahydrofuran- 
naphtalene solution Liu et al. may have extracted not 
only the free lithium, but also the lithium from the LiB y 
compound pushing the concentration of boron towards 
the higher limit of the stability region (y m in,ymax), and 
maybe beyond it. 

To help determine the composition of the LiBj, com- 
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FIG. 7: (color online). Calculated optimized lattice parame- 
ters in phases with linear chains of boron. Top panel: c-axis 
distance (the doubled Li-Li interlayer spacing, 2cLi-u) as a 
function of the B to Li ratio in LiB H . Bottom panel: the 
black circles are the B-B bond lengths in the LiB y phases 
(ffls-s = c/y/2); the red squares are B-B bond lengths in the 
B _1//h phases explained in the text. 



pounds from experimental data we plot the fully relaxed 
lattice parameters in Fig. 7. We observe that the c-axis 
undergoes an almost linear expansion with the increase of 
the boron to lithium ratio: c = 0.365 + 2.746?/. Because 
of the 1~2% systematic errors in the bond length calcu- 
lations within the GGA these results cannot be used to 
pinpoint the absolute value of y. However, the variation 
of the lattice parameters as a function of y is expected to 
be much more accurate and allows one to estimate the 
range of concentrations for the experimentally observed 
compounds. For example, the measured c = 2.875(2) 
A and c = 2.792(1) A values corresponded to nominal 
compositions y — 1.0 and y = 0.82, respectively 0|; the 
slope from Fig. 7 indicates that the difference in y in 
the synthesized compounds was, in fact, about 6 times 
smaller (Ay = 0.03). The discrepancy in the measured 
(2.796 A, Ref. 0) and calculated (3.102 A, Ref. 0) 
c-axis values pointed out by Rosner and Pickett [39j| can 
be explained as that the synthesized compound was not a 
stoichiometric lithium monoboride but rather LiBj,~o.89- 
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FIG. 8: (color online). Calculated total density of states 
(DOS) for the Li2«B m compounds (solid black lines) and hy- 
pothetical B _1,/h phases explained in the text (dotted red 
lines). Fermi level is at eV. 
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FIG. 9: (color online). Calculated partial density of states 
(PDOS) for the Li2nB m compounds. The curves correspond 
to B-s (green), B-p 2 (red), B-p xy (blue), and Li (black) char- 
acters. Fermi level is at eV. 



It should be noted that the determination of the con- 
centration (y = o,b-b I CLi-LU c-Li-Li = c/2) presents dif- 
ficulty only because the length of the boron-boron bond 
is hard to extract from the experiment |46| . However, it 
is the covalent boron-boron bond, being very rigid, that 
determines the c-axis dimension. This can be best il- 
lustrated by simulating linear chains of boron in hypo- 
thetical structures B -1 ^ where lithium is replaced with 
the equivalent number of extra electrons 1/y (a uniform 
positive background is used here to impose charge neu- 
trality). The B -1 /^ structures with one boron atom per 
unit cell keep the basal lattice vectors of the correspond- 
ing LiB y compounds while the c-axis is optimized. Fig. 
7 shows that the charge transfer from lithium to boron 
alone can explain the variation of the bond length for 
y = 0.8 ~ 1.0. The origin of the nearly constant off- 
set between the two curves becomes evident when the 
lithium sublattice is, in turn, simulated without boron. 
A hypothetical Li +1 structure with the same in-plane di- 



mensions as in cv-LiB has a much shorter equilibrium o 
axis of 2.22 A. The E Ll +i (c) dependence in the 2.5~3.1 A 
range of c is almost linear with a slope of 0.76 eV/A. As 
a result, the lithium sublattice in the LiB y phases exerts 
a small stress on the linear chains of boron for all y from 
1 to 0.8. The stress induces the shortening of the boron- 
boron bond by 0.029 A (y = 1) and 0.033 A (y = 0.8) 
and thus turns out to be the main reason for the 0.035 A 
(y = 1) and 0.032 A (y = 0.8) bond length mismatches 
in the simulated phases with and without lithium (Fig. 
7). The result illustrates why LiB y can be represented 
well as a superposition of the two electronically indepen- 
dent doped sublattices. However, our next test shows 
that lithium does play an important role in defining the 
optimal composition of LiB y by affecting the electronic 
states of boron near the Fermi level. 

To further investigate the mechanism leading to the 
existence of the off-stoichiometry lithium borides we plot 
the total and partial DOS for several Li2„B m phases in 
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Figs. 8 and 9 and the band structure in a representative 
a-LiB phase in Fig. 10. The states near the Fermi level 
are hybridized pw-B and Li states 39]. The average pres- 
ence of the Li character in the DOS in the -10 to 3 eV 
energy range is small in a-LiB (N Ll /N^ « 0.3), but be- 
comes more significant in L14B3 (N Ll /N^ w 0.5). The 
pir boron states extend into the lithium-filled interstitials 
the furthest, so they are affected by the electrostatic po- 
tential from the lithium ions the most (in fact, the van 
Hove singularity in a-LiB at E = is not present in the 
corresponding B _1 structure described above, see dotted 
curves in Figs. 8 and 10). Note that the bonding and 
antibonding boron states are not well separated in the 
LiBj, phases: in the case of a-LiB both types are present 
in the 0-1.3 eV energy range (Fig. 10). Since the Fermi 
level in a-LiB (Fig. 8) already catches the edge of the 
antibonding pn-H states, one would naturally expect for 
the compound to benefit from losing lithium. However, 
the system does not follow the rigid band scenario as the 
concentration of lithium increases: the van Hove singu- 
larity is pushed away from the Fermi level (in a-LiB) to 
the right by over 2 eV (in Li 4 B 3 ) (see Fig. 8). The op- 
timal position of the Fermi level near the bottom of the 
pir pseudogap is achieved in LiioBg and LigB7, which 
is a part of the reason why the LiBj, phases have the 
minimum formation energy at y re 0.9. The correlation 
between the position of the Fermi level in the pseudogap 
and the maximum stability has been observed in various 
systems jlli^. 

Comparison of the MS phases with the known metal 
borides phases under pressure. To check whether there 
are more optimal charge transfer and strain conditions 
than those in MSI and MS2 we simulate MS phases at 
other concentrations (Fig. 6, top panel). In the Li-B sys- 
tem the MS3 and MS4 phases have energies well above 
the a-B<-^a-LiB and a-LiB<-^fcc-Li tie-lines, so they will 
be unstable against phase separation into the known com- 
pounds (these MS phases remain metastable under hy- 
drostatic pressure as well) . This test confirms our earlier 
finding that MSl-LiB and MS2-LiB are particularly sta- 
ble because of the near-optimal occupation of the binding 
boron states 0. 

Our conjecture that the hydrostatic pressure would 
favor the MS lithium monoboride phases over the off- 
stoichiometric LiBj, compounds 9j is also supported by 
the results shown in Fig. 6. As far as the possibility of the 
MS-LiB formation is concerned, it is rather unfortunate 
that the LiBj, phases additionally stabilize by becoming 
more lithium rich under pressure. The information in 
Fig. 6 can be used to evaluate the minimum pressure 
required to stabilize MS2-LiB with respect to a-B and 
LiBj,. We estimate by interpolation that at P m in ~ 5 
GPa the MS2-LiB phase lies on the line that goes through 
x = and is tangent to the LiBj, parabola. Therefore, 
MS2-LiB could appear in the experiment at P > P m in 
(not accounting for the possible systematic errors and 
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FIG. 10: (color online). Band structures of a-LiB f|3^. 
solid black lines) and hypothetical B _1 (described in the text, 
dotted red lines) phases. To facilitate the comparison of the 
band structures in the two phases we doubled the B _1 unit 
cell along the c-axis and used the equilibrium lattice constants 
of a-LiB. With the optimized c-axis the B _1 electronic states 
in the -11 ~ -2 eV range upshift by about 0.15 eV. Fermi level 
is at eV. 



the thermodynamics effects |9j) if all other phases in the 
Li-B system were metastable under such pressures. Us- 
ing the parabola coefficients for the three pressures given 
in Fig. 6 we also find by interpolation that MS2-LiB has 
the same formation enthalpy as the most stable LiBj,~o.82 
phase at about 12 GPa. Finally, we observe that the line 
connecting MS2-LiB and fcc-Li crosses the parabolas in 
all the cases for pressures below 30 GPa 50] . This im- 
plies that if the MS2-LiB phase was synthesized, LiB y 
might still be present in the sample as a by-product. For 
analysis of the Li-B system at higher pressures one needs 
to take into account that pure lithium undergoes phase 
transformations from fee to hRl and eventually to c/16 
near 40 GPajEJ. 

The chances for the formation of the MS-LiB phases 
depend on where they are located on the phase dia- 
gram relative to LiBj, and the most lithium rich stable 
phase below 50% concentration. The known phases in 
this region have small atomic volume {Vl^Bu = 6.4 and 
VliB 3 = 7.5 A 3 /atom under ambient conditions |l3[ IsfiL 
H5) compared to MS2-L1B (11.2 and 6.7 A 3 /atom at 
and 30 GPa, respectively). The boron-rich phases could 
potentially bar the formation of MS-LiB under pressure, 
however they would need to have a very low formation 
enthalpy; for example, at P = 12 GPa it would need to be 
below H M S2-LiB = H LiBysl0 g2 = -0.38 eV/atom. Simu- 
lation of a known phase with CaBg prototype, present 
in the K-B system^]], El- could give information on 
how boron-rich phases respond to hydrostatic pressure. 
However, the ordered CaBg-LiBg has a large atomic vol- 
ume and actually becomes less stable under pressure: 
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FIG. 11: (color online). Calculated phase diagrams for Na-B 
and Rb-B systems at different pressures. Hollow red points 
are M2 n B m phases, from left to right: M2B4, M2B3, M4B5, 
M 6 B 7 , q-MB, /3-MB, and M g B 7 . The boron rich phases are 
Na3B2o and CaBs prototypes with Pearson symbols 0.546 and 
cP7, respectively. The vertical dotted lines mark the 50% 
composition. 
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FIG. 12: (color online). The same as in Fig. 11 for K-B and 
Cs-B systems. 



6. OTHER MONOVALENT AND SOME 
HIGHER- VALENT METAL BORIDES 



Hf = 0.016 eV/atom, V = 10.1 A 3 /atom at zero pres- 
sure and Hf = 0.303 eV/atom, V — 8.6 A 3 /atom at 
30 GPa pressure. Rosner and Pickett pointed out that 
a compact pseudodiamond phase B32 (NaTl prototype) 
might appear under pressure 39] . Using enthalpy versus 
pressure curves for the a-LiB and B32 phases we find the 
crossover pressure to be 22 GPa. MS2-LiB stays below 
B32-LiB until about 65 GPa. Overall, our simulations 
suggest that there might be a window of pressures at 
which the MS-LiB phases can be synthesized. 



Alkali and transition metal borides. According to the 
experimental databases and the latest review of alkali- 
metal borides [HEl 

there are no stable sodium or potas- 
sium borides above 15% concentration of metal and no 
stable rubidium or cesium borides in the whole concen- 
tration range. Na3B2o (Pearson symbol 0546) and KBg 
(Pearson symbol cPl) compounds, made out of boron 
polyhedra intercalated with alkali atoms, are considered 
the most metal-rich known borides in the Na-B and K- 
B systems, respectively 0, Our simula- 

tions confirm that these compounds have negative for- 
mation enthalpies of -58 meV/atom for Na3B2o and -29 
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meV/atom for KBq with respect to a-B and fcc-M 35] . 
The formation Gibbs free energies for these compounds 
might be less negative if they were evaluated with re- 
spect to (3-B at finite temperature j^. This could be the 
reason why there is no conclusive evidence of potassium 
hexaboride synthesis 01- Complete theoretical investiga- 
tion of the boron-rich compounds such as Na3B29 (Ref. 
|ft4j ) is beyond the scope of this study but it is interesting 
to see where the nearly stoichiometric MB y phases place 
with respect to the known phases in these systems. 

As in the LiB y compounds, the metal and boron sub- 
lattices in the alkali boride phases MB y (M = Na, K, Rb, 
Cs) are found to be very weakly correlated. In fact, the 
larger alkali atoms push the boron chains farther apart 
(see Table I) weakening the boron interchain bonds (note 
the lower energy difference between q-MB and (3-MB 
listed in Table I) . This again leads to the situation when 
MB y (M = Na, K, Rb, Cs) compounds can easily adapt to 
an optimal composition by having incommensurate metal 
and boron sublattices. Figures 11 and 12 demonstrate 
that in all Na-B, K-B, Rb-B, Cs-B systems the stoichio- 
metric phases with linear chains of boron prefer to lose 
some metal, the opposite tendency compared to the Li-B 
system. While the most stable LiB y composition appears 
to be determined primarily by the optimal level of boron 
doping, in the larger alkali- metal borides, MB y , the lat- 
tice mismatch between the metal and boron sublattices 
must be playing a more significant role. Note that the 
formation enthalpy points for these Ma n B m phases are 
not symmetric, curving up more rapidly in the metal rich 
region. 

The stabilization from losing a few percent of alkali 
metal is noticeable but not enough for NaB y and KB y 
to have a negative formation enthalpy at zero pressure. 
Because the interchain spacing is determined mostly by 
the alkali cations the C\\ + C12 force constant in (3-MB 
decreases as one moves down the periodic table (see Ta- 
ble I) . C33 also becomes smaller as the boron-boron bond 
length gets longer. The softness of the MBj, phases (M = 
Na, K, Rb, Cs) invites the use of hydrostatic pressure for 
their synthesis. Moreover, in the Na-B and K-B systems 
the MB y phases stabilize more rapidly than the MS2-MB 
phases, which makes it unlikely for the latter to form un- 
der the pressures considered. Synthesis of MB y (M = 
Na, K, Rb, Cs) would provide valuable information on 
the ways the linear chains of boron could be stabilized. 
Because the alkali-metal borides are not fully explored, it 
would not be surprising if a not considered here or a com- 
pletely unknown phase appeared in such an experiment. 

While the /3-MB phases are only metastable for the 
borides in the alkali-metal series, there are two reported 
stable transition-metal monoborides in this configura- 
tion: RhB and PtB (Ref. [HEi]). Our fully relaxed unit 
cell parameters (see Table I) agree well with experiment 
for /3-RhB (a = 3.309 A, c = 4.224 A), but they dis- 



agree by over 10% with the measured values (a = 3.358 
A, c = 4.058 A) for /3-PtB [H HU . Identification of the 
source of this discrepancy requires additional study of 
this system. Nevertheless, the data on /3-RhB and (3- 
PtB in Table I give an idea about what difference the d- 
electrons cause in the boron-boron binding compared to 
the case of the alkali-metal monoborides. For example, 
/3-RhB and /3-PtB no longer have the optimally doped 
double-bonded boron chains: the boron-boron bond is 
so overstretched that it exceeds the sp 2 bond length in 
the AlB2-type compounds, resulting in the increase of 
the c-axis compressibility compared to the alkali-metal 
monoborides. The significant reduction of the interchain 
distances leads to an over 300% increase in the C\\ + C12 
force constants. One more important consequence of the 
more compact arrangement of atoms in the lateral di- 
rection and the hybridization of the <i-orbitals of metal 
with the valence states of boron is the much larger en- 
ergy difference between the (3-MB and a-MB structures. 
This makes the formation of the off-stoichiometry phases 
with linear chains of boron in the transition-metal mono- 
borides energetically unfavorable. 

Noble- and divalent-metal borides. AgB2 and AuB2 
have been shown to have big positive formation 
energies 4J , so they should not form at ambient pressure. 
Recent experiments suggest that some superconduc ting 
Ag-B phase was formed by pulsed laser deposition[l2|. 
Because of the synthesis conditions, the thin-film sam- 
ples were inhomogeneous and did not produce new x-ray 
peaks. It was assumed that the observed phase was |a|- 
prototype, although the T c turned out to be much lower 
than the anticipated value^l). Formation of Ag vacan- 
cies in |a|-Aga;B2 was suggested by Shein et al. as a pos- 
sible explanation of the observed data^|- Current sim- 
ulations offer other possibilities for formation of phases 
under non-equilibrium conditions: the proposed phases 
still have positive formation energies but they are less 
unstable and below the respective \a\ <-> metal tie-lines 
(Fig. 13). The intermediate phases in the concentration 
range from 33% to 60% generally stay below the tie- lines. 



TABLE I: Calculated properties of /3-MB metal borides 
(NiAs prototype). 



compound 


a 


CO 


Cn + C12 


C33 


Efi-MB — E a -MB 




(A) 


(A) 


(GPa) 


(GPa) 


(meV/atom) 


LiB 


4.013 


3.120 


139 


542 


-10 


NaB 


4.697 


3.196 


111 


379 


-3.1 


KB 


5.390 


3.240 


79 


268 


-1.6 


RbB 


5.662 


3.267 


77 


251 


-1.7 


CsB 


5.976 


3.325 


69 


252 


-2.0 


RhB 


3.382 


4.185 


498 


309 


-267 


PtB 


3.765 


3.655 


546 


293 


-108 
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FIG. 13: (color online). Calculated phase diagrams for mag- 
nesium and noble-metal borides. The solid points are layered 
phases from left to right: \a\, MS4, MS2, MS3. We also show 
5-AuB2 phase ISl^l. which is significantly lower than |a[-AuB2. 



form and could possibly exist only in the form of a de- 
fect in |a|-MgB 2 . We calculate the following series: \a/3\, 
\a\aft\, \ot\ot\otj3\, and so on up to 8 |a| unit cells, and find 
by extrapolation that the energy required to insert a sin- 
gle magnesium layer into the a|-MgB 2 matrix is quite 
high: over 25 meV per atom in the additional layer of 
magnesium. 



7. SUMMARY 

Electronic properties. It is illustrative to compare the 
important features of the electronic structure in the low- 
valent metal borides to those in the lithium borides, 
which were discussed in our previous study 9]. We focus 
mainly on the p-states of boron, important for the stabil- 
ity and the superconductivity in these compounds. For 
convenience, we calculate the band structure and PDOS 
for \a\, MS2, and MS3 phases because all three have a 
hexagonal unit cell. The key characteristics of the boron 
states, along with parameters of the unit cells for these 



If formed, MS2-AgB would be even more dynamically 
stable than MS2-LiB: we find that the former represents 
a stable equilibrium with the frequencies of the softest 
interlayer modes u/ x , y three times higher than those in 
the latter 0. This is an expected result, considering that 
the silver bilayer in MS2-AgB remains bound by the d- 
electrons even if it donates most of the charge from the 
s-orbital to boron as lithium does. The electronic prop- 
erties and the stability of the MS noble-metal borides are 
further discussed in Section 7. 

We find that the o-MB 2 phasej^ is surprisingly much 
more stable than the |cv|-MB 2 phase for several metals: 
Au, Ca, K, Pd, and Pt (by 215, 28, 214, 28, and 272 
meV/atom, respectively). Apparently, metal atoms pre- 
fer to hybridize their valence states with p-orbitals of 
boron more strongly by shifting to a boron-boron bond, 
rather than simply donate their valence electrons. Note 
that the five metals are either large in size (dc a ^3. 88 
A, dff=4.71 A, see Fig. 4) or have a big work func- 
tion (bulk values: 0^ u =5.1 eV, (f>pd=5.12 eV, </>p t =5.65 
eV|5o. The discovery of the lower-symmetry 5- 

AuB 2 phase |32]| rules out the possibility of cv|-AuB 2 syn- 
thesis. The formation energy of <5-AuB 2 remains positive 
(0.45 eV/atom, see Fig. 13), which makes this phase un- 
likely to form as well. In the Ca-B system <5-CaB 2 is still 
unstable against phase separation into GslBq and fcc-Ca 
0E3 b y 144 meV/atom. 

The MS phases in the Mg-B system stay at least a few 
meV/atom above the |a|-MgB 2 ^->hcp-Mg tie-line. Hy- 
drostatic pressure is insignificant to their relative stabil- 
ity because they compete against similar phases. There- 
fore, the MS magnesium boride phases are not likely to 



TABLE II: Metal boride phases: formation energy (Ef, 
eV/atom), relative stability (AEx, eV/atom, Eq. 1), in- 
plane boron-boron bond (clb-b, A), interplanar boron-boron 
distance (cb-b, A), internal coordinate of the metal atom 
(zm ^l^ftl. l2fll |. position of pa-band in boron at T and A k- 
points (E r , E A , eV), and PDOS at E F for B-p xy (N^ y (0), 
states/(eV-spin- boron atom)). 



phase 


Li 


Ag 


Au 


Cu 


Mg 




Ef 


0.003 


0.520 


0.667 


0.371 


-0.131 




CLB-B 


1.717 


1.744 


1.737 


1.721 


1.776 


M 


CB-B 


3.469 


4.081 


4.260 


3.382 


3.521 


(33%) 


Er 


1.48 


1.29 


1.26 


1.14 


0.39 




Ea 


1.70 


1.01 


0.95 


0.65 


0.77 




a£„(°) 


0.076 


0.090 


0.093 


0.106 


0.049 




Ef 


-0.162 


0.374 


0.488 


0.266 


-0.089 




AEmS2 


-0.164 


-0.016 


-0.012 


-0.013 


0.009 


MS2 


&B-B 


1.765 


1.739 


1.734 


1.734 


1.805 


(50%) 


CB-B 


5.522 


6.369 


6.589 


4.915 


5.989 


m 


ZM 


0.496 


0.368 


0.356 


0.388 


0.440 




E T 


0.99 


1.19 


1.17 


0.69 


0.15 




Ea 


0.99 


1.19 


1.17 


0.69 


0.15 




K v (o) 


0.059 


0.086 


0.091 


0.091 


0.043 




Ef 


-0.117 


0.305 


0.394 


0.277 


-0.061 




AEmss 


-0.118 


-0.007 


-0.006 


0.054 


0.018 


MS3 


clb-b 


1.745 


1.731 


1.722 


1.695 


1.816 


(60%) 


CB-B 


8.318 


8.785 


9.097 


6.903 


8.511 


f29] 


ZM 


0.178 


0.231 


0.237 


0.228 


0.194 




Er 


1.23 


1.22 


1.25 


0.88 


0.03 




Ea 


1.23 


1.22 


1.25 


0.89 


0.03 




a£„(o) 


0.066 


0.090 


0.098 


0.074 


0.039 
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phases are given in Table II. 

For all MS2 metal borides the pa band along T-A is 
practically flat because of the large separation between 
boron layers, as shown in Fig. 4 of Ref.Q for MS2-LiB 
and in Fig. 14 for MS2-AgB. For Ag and Au this band 
does not move much from the respective average posi- 
tions in \a\ and the PDOS of pa states in boron at the 
Fermi level, (0) , stays nearly the same (Table II and 
Fig. 15). In fact, these boron properties remain the 
same in the more metal-rich structure MS3. The results 
suggest that the level of doping of the boron layers in 
the MS noble metal borides is nearly independent of the 
number of metal layers. Considering this and the fact 
that the MS phases in the 33%-100% range closely follow 
the | a | <->fcc lines on the Ag-B and Au-B phase diagrams 
these MS phases can be viewed as a mixture of weakly 
interacting building blocks: the |a| unit cell with an es- 
tablished charge redistribution within it and the closed- 
packed layers of pure metal. Therefore, one could ex- 
pect the superconducting properties of the boron layer 
in the hypothetical phases |a|-AgB2 and the MS silver 
borides in the 33%-100% concentration range to be sim- 
ilar. However, the nonequilibrium conditions necessary 
to synthesize such compounds^3| may introduce disorder 
destroying the superconducting states Q- 

Compared to the MS gold and silver borides, the MS 
copper borides deviate more from the |a|-CuB2<->fcc-Cu 
line (Fig. 13) and the level of boron doping is more non- 
monotonic as a function of the metal concentration (Ta- 
ble II) , even though the work functions of Cu and Ag are 
close (bulk values: (f>cu — 4.65 eV, <f>A g — 4.26 eV, Ref. 
[5(|). The different behavior must be the result of the 
more pronounced strain between the copper and boron 
layers caused by the smaller size of copper (dcu = 2.563 
A, d Ag = 2.937 A, d Au = 2.949 A, see Fig. 4). In the 
more electron-rich Mg-B system the level of boron dop- 
ing increases as additional magnesium layers are added. 
In MS3-Mg3B2 the pa states are almost completely oc- 
cupied. 

Lithium turns out to be a special case among the mono- 
valent metals: it has the right size and can easily do- 
nate electrons to stabilize the boron layers in the MS 
phases. Fig. 6 and Table II show that the stabilization 
effect is largest for the 50% concentration MSI and MS2 
lithium borides. No less importantly, by giving up most 
of their charge the lithium layers interact only weakly 
with one another, which makes the compound very soft 
along the c-axis and gives the opportunity to stabilize the 
compound even more with hydrostatic pressure [(|. The 
addition of extra layers does not result in further sta- 
bilization, as boron in MS3-Li3B2 is doped less than in 
MS2-LiB, judging by the location of the pa band along 
r-A (Table II and Fig. 15). The present PAW calcu- 
lations confirm our previous result that MS2-L1B has a 
higher PDOS of boron pa states at the Fermi level than 
that in MgB 2 (Ref. 0). Note that the PDOS in our sim- 
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FIG. 14: (color online). Band structure and partial density 
of states (PDOS) in a hy pothetical MS2-AgB phase, calcu- 
lated in APW+lolU Hj. PDOS units are states/(eV-spin) 
per boron atom. The thickness of band structure lines is pro- 
portional to boron p x<y (red) andp z (green) characters. Fermi 
level is at eV. 
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(green) states in lithium and silver MS borides calculated with 
VASP. Fermi level is at eV. 



ulations is found by decomposition of the wavefunction 
within a sphere of fixed radius and can slightly vary with 
this parameter, as well as with the approximation used. 
In the APW+lo calculation we obtained a 12% increase 
of Nj? (0) in MS2-L1B compared to that in MgB 2 using 

Rmt = 1-6 a - u - Mi while in the present PAW calculations 
we observe a 20% increase using the default PAW radius 
of 1.7 a.u. 
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The ample amount of the boron pa states at the Fermi 
level in MS2-LiB holds great promise for this compound 
to be a good superconductor. However, a more thorough 
calculation of the electron-phonon coupling, such as in 
Ref. is required to say with certainty whether the 

new lithium monoboride can compete with the record- 
holder MgB 2 . Such calculation is under way |6l). 

X-ray reflections. As we pointed out previously 0, the 
resulting structure in metal borides at 1:1 composition 
could be a random mixture of MSI and MS2, because 
they differ only by a long-period shift in stacking order 
and, therefore, are nearly degenerate. The two phases 
also have very close nearest neighbor distances in most 
systems and random structures would still have a con- 
stant separation between boron layers. For magnesium 
and silver MS monoborides these periods correspond to 
29 = 16.4° and 26 = 16.8°, respectively (for A = 1.5418 

A). 

The largest difference (3%) in the c-axis for MSI and 
MS2 is actually found for Li-B, resulting in 29 = 16.6° 
{dusi-LiB = 5.35 A) and 29 = 16.1° {&ms2-ub = 5.52 
A); the more lithium-rich phase MS3-Li3B2 would pro- 
duce a reflection at 29 = 10.6° {dMS3-Li 3 B 2 — 8.32 
A). The most pronounced peaks in the published x-ray 
data for lithium monoboride are at 29 = 25.5°, 41.3°, 
and 45.0° which fit well to the calculated a-LiB x-ray 
pattern]^. Interestingly, two reflections at low angles 
29 = 12.2° and 20.9° were observed at 40-50% of lithium 
concentration^]- A low-angle reflection was also de- 
tected at 12.8° in Ali_ a; Li a ;B2 under heavy Li doping^. 
However, none of the observed peaks in the samples pre- 
pared at ambient pressure match the calculated x-ray 
reflections in the MS lithium borides. 



8. CONCLUSIONS 



composition. The MS lithium monoboride phases are 
shown to have lower formation enthalpy with respect to 
the experimentally observed nearly stoichiometric LiBj, 
phases under hydrostatic pressure. This encouraging 
result suggests that the new superconducting MS-LiB 
phases might form under proper conditions. The lowest 
required pressure depends on the position of boron-rich 
phases in the Li-B phase diagram and could be as low as 
several GPa. 

iv. For a more complete description of the Li-B sys- 
tem we introduce a simple model of the off-stoichiometric 
LiBj, phases which explains the available experimental 
data. We demonstrate that because of the weak cor- 
relation between the boron and lithium sublattices the 
compound can easily adapt to an optimal composition, 
which corresponds to an optimal level of boron doping 
with the Fermi level lying near the bottom of the pseu- 
dogap. Interestingly, these phases turn out to be stable 
over a range of concentrations around y = 0.9, in excel- 
lent agreement with experiment. We list the relaxed unit 
cell parameters which should be helpful in the determi- 
nation of the LiBj, composition. 

v. We consider the MB y phases for other alkali-metal 
borides and find that these compounds also benefit from 
going off-stoichiometry, only in this case they prefer to 
lose some metal. The ensuing gain in enthalpy does 
not make them stable under ambient conditions, however 
these phases might form under hydrostatic pressure. Syn- 
thesis of the MBj, phases (M = Na, K, Rb, Cs) would pro- 
vide valuable information on how linear chains of boron 
could be stabilized. 

We thank M. Calandra, F.H. Cocks, V. Crespi, P. Lam- 
mert, E. Margine, and J. Sofo for valuable discussions. 
We acknowledge the San Diego Supercomputer Center 
for computational resources. 



The main results of the present study can be summa- 
rized as follows: 

i. We have identified a previously unknown class of 
metal-rich layered phases that are comparable in energy 
to existing metal borides. This interesting accidental re- 
sult should be credited to the exhaustive consideration 
of all candidates in the DMQC method and the careful 
structural relaxation in the calculation of their ground 
state energies. 

ii. Our ab initio results suggest that the MS phases are 
most suitable for electron-deficient metal boride systems. 
In the Ag-B and Au-B systems the MS phases are less 
unstable than the corresponding diborides with the A1B 2 
prototype but they still have positive formation energies. 
In the Mg-B system the MS phases are metastable and 
could possibly exist only as a defect in MgB2. 

iii. The MS-LiB phases present a special case among 
the MS metal borides: lithium has the right size and 
valence to stabilize the hexagonal layers of boron at 1:1 
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